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Abstract 

^ . Comparing proposed quantities to analyze charged particle fluctuations in heavy ion ex- 

periments we flnd the dispersion of the charges in a central rapidity box as best suited. Various 
energies and difi^erent nuclear sizes were considered in an explicit Dual-Parton-Model calcu- 
lation using the DPMJET code. A definite deviation from predictions of recently considered 



(N 

en 

^— «, ' statistical models was obtained. Hence the charged particle fiuctuations should provide a 

f~^ , clear signal of the dynamics of heavy ion processes. They should allow to directly measure 

f^ ■ the degree of thermalization in a quantitative way. 

^' 1 Introduction 

O . In the analysis of the hadronic multi-particle production (for a recent review see ||l|) a key ob- 

servation has been the local compensation of charge|g). The charge fluctuations connected to the 
soft hadronic part of the reactions were found to involve only a restricted rapidity range. This 
observation limited the applicability of statistical models to rather local fluctuations (see e.g. ||^ 
where the charge fluctuation between the forward and backward hemisphere^] was discussed). 

C^ ' In heavy ion scattering charge flow measurements should be analogously decisive. It is a central 

question of an unbiased analysis whether the charges are distributed just randomly or whether 
there is some of the dynamics left influencing the flow of quantum numbers. This is not an 
impractical conjecture. In heavy ion experiments the charge distribution of particle contained in 
a central box with a given rapidity range [— l^nax. , +Fmax.] can be measured and the dispersion of 
this distribution: 

< SQ^ >=< [Q- <Q>f> (1) 

can be obtained to sufficient accuracy. For sufficiently large gaps this quantity also reflects the 
long range charge flow. 

It was proposed to use this quantity to distinguish between particles emerging from a equi- 
librized quark-gluon gas or from a equilibrized hadron gas[|[ |[ 0- In a hadron gas each particle 
species in the box is taken essentially poissonian. In a central region where the average charge 
flow can be ignored, one obtains a simple relation for particles like pions with charges and ±1 : 

< SQ^ > = < TVchargod > . (2) 

It is argued in the cited papers that this relation would change in a quark gluon gas to: 

< SQ^ >= J2 9^' < ^' >= 0-19 < ^charged > (3) 



where qi are the charges of the various quark species and where again a central region is consid- 
ered. The coefficient on the right was calculated|g] for a two flavor plasma in a thermodynamical 
consideration which predicts various quark and gluon contributions with suitable assumptions. A 
largely empirical final charged multiplicity iVchaigcd = |(^giuc + l-ZA'^quark + l-2A^antiquark) was 
used. 

It should be pointed out that the estimate is not trivial. In the considered mquark = theory 
the observable Q is not infrared safe and makes no sense. A way to make it well defined[H is to 
consider the quantity: _ 

^Jquark ^ ^iquark ^::quai'k ^ ^/quark ^ \^) 

which avoids the infiuence of extra sea quarks. Such a correction is also needed to have an 
observable which can be expected to survive hadronization. Numerically the effect is not very big 
and the problems can be ignored if only a rough description is sufficient. 

There are a number of sources of systematic errors in the above comparison. The result 
strongly depends on what one chooses as primordial and secondary particles. Considering these 
uncertainties we follow Fialkowski's conclusion ||g] that a clear cut distinction between the hadron- 
and the quark gluon gas is rather unlikely. This does not eliminate the interest in the dispersion. 
The hadron gas model is anyhow no optimal reference point to compare with. 

In the next section we discuss various possible measures to observe such fluctuations. We favor 
the dispersion of the charge transfer. Using an expHcit Dual Parton model calculation we observe 
a clear distinction between models with local compensation of charge and equihbrium approaches. 
In section 3 a simple interpretation of the dispersion in terms of quark lines is outlined. This 
suggests to compare the dispersion to the particle density as it is done in section 4. 

2 Various Measures for Fluctuations in the Charge Distribu- 
tion 

For the analysis of the charge structure several quantities were discussed in the recent literature. 
It was proposed to look at the particles within a suitable kinematic region and to measure just 
the mean standard deviation of the ratio R of positive to negative particles: 

N+ N. ^^ 



or the quantity F: 



<<5i?^>=<(_r_<_r>) > (5) 
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<SF'>^<(-^ <^7^—>] > (6) 

\ -^* charged J* charged / 

where Q = N+- N-. 

We consider them not attractive. The quantities are not suitable for small intervals, as there 
are actually undefined in a certain region. They are less clean than the dispersion, < SQ >, as 
they are not exclusively connected to the fiavor structure and as they mix up charge and density 
fiuctuations. 

For large nuclei at high energies this is not a problem as the density fiuctuations are small and 
all three quantities are connected by the following relations JG): 

< A^chargcd >< 5R^ >= 4 < A^chargcd >< 5F^ >= 4 • ^^ ^ (7) 

^ -^^chargcd ^ 

To examine the range where these relations hold, all three quantities were calculated in the Dual 
Parton model implementation DPMJETpH] . For the most central 5% Pb-Pb scattering at LHC 
energies {^/s = 6000 A GeV) there is indeed a perfect agreement between all three quantities as 
shown in figure 1. For the most central 5% S-S scattering at RHIC energies {^/s — 200 A GeV) 
the agreement is no longer as good. For the minimum bias S-S scattering at RHIC energies the 
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<Nch>< 5R^ > LHC 

4 < N,h >< SF^ > LHC 

4<5Q^> / <N,h> LHC 

i<5(^ > / <Nah> RHIC 

= 4 < 5Q2 > / < 5V,A > SPS 




Figure 1: Charge fluctuations for the most central 5% Pb-Pb scattering at LHC energies (^/s = 
6000 A GeV). Also shown are dispersions for RHIC (^/i = 200 A GeV) and SPS (Eab A GeV) 
energies. 

agreement is lost and the first two expressions behave rather erratic. The same erratic behavior is 
seen for the proton-proton case which is shown figure 2. As any conclusion will strongly depend on 
a comparison of central processes with minimum bias events, we advocate to stick to the dispersion 
of the net charge distribution < SQ"^ >. 

We observed no significant difference between rapidity and pseudo-rapidity boxes. 



3 A Simple Relation between the Quark Line Structure and 
Fluctuations in the Charge Flow 

To visualize the meaning of charge fiow measurements it is helpful to introduce a general factor- 
ization hypothesis. It is not exact but it is widely expected to hold to good accuracy. It postulates 
that the flavor structure of an arbitrary amplitude can be described simply by an overall factor, 
in which the contribution from individual quark lines factorizeQ 

The hypothesis can be used to obtain the following generalization of the Quigg-Thomas relation p^, 
[l3| , Q. This generaHzation states that the correlation of the charges Q{yi) and Q(y2), which are 

1 The hypothesis is based on the exchange degeneracy of octet and singlet Regge trajectories effectively changing 
the 5f7(Af(,ffoctivo) flavor symmetry to an (/(Afcffoctivo) symmetry in which this relation is exactly valid. One of the 
important corrections to the hypothesis originates in the special behaviour of the masses of the lowest lying mesons 
of the trajectories, which is especially significant in the pseudo-scalar sector between the ttq and the rj meson. 
That one of the two neutral states is sometimes suppressed, introduces a mild anticorrelation between neigbouring 
flavors, which can be ignored for our consideration centered at long range charge transfers. 

If a higher accuracy is desired the hypothesis can be restricted to primordial particles generated at a high 
"temperature" and "secondary" charges produced during the decay of primordial particles can be considered extra. 
With simple assumptions their contribution can be related to the corresponding particle spectrum. If all charged 
particles were secondaries the dispersion of the charge transfer across an arbitrary rapidity boundary would be given 
by the Quigg-Thomas relation|l(j, h]|, rLa| < SQ^(y) >= o"^P(,harged(y) where a = 1 if widening and narrowing effects 
balance. 
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Cq=< nL X^R' > 10 TeV -•- 
CQ = i<Kh>< SF^ > 10 TeV --A- 
Cq = 4 < (5g2 > / < N,h > 10 TeV --^- 



Co = 4 < jg^ > / < Af,6 > 200 GeV 
Cq^4<Iq^ > / < N^h > 20 GeV 




Figure 2: Charge fluctuations for minimum bias pp scattering at SPS, RHIC and LHC energies 



exchanged during the scattering process across two separate kinematic boundaries, is just: 

< {Q(2/l)~ < Qiyi) >} {(9(2/2)- < Q{y2) >} >= ncommonlincs < (q- < Q >f > ■ (8) 

where ^common lines counts the number of quark lines intersecting both borders and q is the charge 
of the quark on such a hne. Depending on the flavor distribution average values < {q— < q >)^ = 
0.22 •••0.25 are obtained. 

Most observables of charge fluctuations depend on this basic correlation. The fluctuation of 
the charges within a [— ^max., +i^nax.] box discussed above contains a combination of three such 
correlations: 



<SQ^>^< SQiyif > + < 5Q{y2f > - 2 < SQ{yi) ■ SQ{y2) > 
Using (|8|) the dispersion of the charges in a box subtracts to: 

< SQlhoxf >= niinos entering box < (g- < 9 >)^ 

where nunes entering box is the number of quark lines entering the box. 



(9) 



(10) 



4 Calculation of the Dispersion of the Charge Distribution 
within a Box 

Let us consider the prediction of a thermodynamic model in more detail. In the thermodynamic 
limit with an infinite reservoir outside and a finite number of quarks inside all quark lines will 
connect to the outside as shown in figure 3. The dispersion of the charge transfer is therefore 
proportional to the total number of particles inside. In the hadron gas all particles contain two 
independent quarks each contributing with roughly 1/4 yielding the estimate of (||). For the 
quark gas only one quark originates in the thermalization and is taken to be responsible for charge 
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Figure 3: Quark lines entering the box in the thermodynamic limit 
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Figure 4: Quark lines entering the box with local compensation of charge 



transfer yielding a considerable smaller result. Obviously there are several refinements to this 
simple picture. 

Let us consider the limit of a tiny box. Looking only at the first order one trivially obtains: 



< (5g2 > / < ATchargcd >= 1 



(11) 



which corresponds to the hadron gas value. 

If the box size increases to one or two units of rapidity on each side this ratio will typically 
decrease, as most models contain a short range component in the charge fluctuations usually 
attributed to secondary interactions. One particular short range fluctuation might be caused by 
the hadronization of partons of the quark gluon gas. The quark antiquark pair needed for the 
hadronization is assumed to be short range so that for a box of a certain size one just obtains the 
charge dispersion of the original partons. This is responsible for the reduction of the fluctuation 
discussed above. The decreasing is however not very distinctive. Common to many models are 
secondary interactions which involve decay processes and comover interaction. In hadron hadron 
scattering processes such correlations are known to play a signiflcant role and there is no reason 
not to expect the same for the heavy ion case. 

After a box size passed the short range the decisive region starts. In all global statistical 
models!^, |^, ^ the ratio will have to reach now a flat value. Only for very large rapidity ranges 
charge conservation will force the ratio to drop (oc 1 — j/max./^in.max.)- This is different in string 
models as it is illustrated in flgures 1 and 2 . The model calculation shows with its rapid fall 
off a manifestly different behavior. It is a direct consequence of the local compensation of charge 
contained in string models. The effect is illustrated in Figure 4 in which only quark Hues which 
contribute to the charge flow and which intersect the boundary are shown. The local compensation 
of charge allows now only a contribution of lines originating around the boundaries. If the distance 
is larger that the range of charge compensation the dispersion will no longer increase with the box 
size. The total contribution will now be just proportional to the density of the particles around 
the boundaries 

< 60^ > ex Pcharged(ymax.)- (12) 

This resulting scaling is illustrated in a comparison between both quantities shown in Figure 5 
for RHIC and LHC energies. The agreement is comparable to the proton proton case shown 
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Figure 5: Comparison of the dispersion of the charge distribution with the density on the boundary 
of the considered box for central gold gold resp. lead lead scattering at RHIC and LHC energies. 

in Figure 6. The proportionality is expected to hold for a gap with roughly ^Sy > 1 . For 
smaller boxes some of the quark Hnes seen in the density do not contribute as they intersect both 
boundaries. For large rapidity sizes there is a minor increase from the leading charge flow Ql 
originating in the incoming particles. In a more careful consideration p^ one can subtract this 
contribution 



<. Oi^ ^leading charge migration ~ 

and concentrate truly on the fluctuation. 
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(13) 



The prediction for the proportionality factor for the case of mere short range fluctuations would 
be roughly a factor one (see footnote 1). In string models primordial particles are responsible 
for a longer range charge transfer coming from the contributions of the quark resp. diquark 
fragmentation chains. Taking everything together one obtains 



< SQ"^ >= 22 ^ ?^strings- < (q- < q >f > + /secondary CT - Pcliargod(y) } 

left+right 



(14) 



where ngtrings is the number of strings where /secondary is the fraction of secondary to primary 
particles and where the width of the local fluctuations a is roughly unity. In our Dual Parton 
Model calculation we observe factor of roughly 1.2 between the density and the dispersion. It 
means that most of the fluctuation originate in secondary interaction and that the effective larger 
coefficient of the flrst term which signiflcantly rises the factor plays only a lesser role. 



Conclusion 

In the paper we demonstrated that the dispersion of the charge distribution in a central box of 
varying extend in rapidity is an extremely powerful measure. Within the string model calculation 
the dispersion seen in relation to the spectra shows no difference between simple proton proton 
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Figure 6: Comparison of the dispersion of the charge distribution with the density on the boundary 
of the considered box for proton proton scattering at RHIC and LHC energies. 

scattering and central lead lead scattering even though both quantities change roughly by a factor 
of 400. 

The dispersion allows to clearly distinguish between conventional string models and thermal 
models. As most expected changes in the dynamic are somehow connected to the onset of thermal- 
ization they will introduce more fluctuations. Even if the truth should He somewhere in between 
string models and quark gluon plasma models it is therefore quite reasonable to hope that the 
position of the transition can be determined in a quantitative way. 
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